Tuning of tunneling current noise spectra singularities by localized states charging 
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We report the results of theoretical investigations of tunneling current noise spectra in a wide 
range of applied bias voltage. Localized states of individual impurity atoms play an important role in 
tunneling current noise formation. It was found that switching "on" and "ofF of Coulomb interaction 
of conduction electrons with two charged localized states results in power law singularity of low- 
frequency tunneling current noise spectrum (!//") and also results on high frequency component 
of tunneling current spectra (singular peaks appear). 
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I. INTRODUCTION 

In the present work we discuss one of the possible rea- 
sons for the noise formation in the wide range of applied 
bias in the STM/STS junctions. Up to now the physical 
nature and the microscopic origin of tunneling current 
noise spectra formation in low and high frequency region 
is in general unknown. Only the limited number of works 
was devoted to the problem of 1//" noise study and we 
have found only a few works where high frequencv region 
of the tunneling current spectra is studied [i| , 0] ■ 

We suggest the theoretical model for tunneling current 
noise spectra above the flat surface and above the impu- 
rity atoms on semiconductor or metallic surfaces. Our 
model gives an opportunity to describe not only singular- 
ities in a low frequency part of tunneling current spectra 
but also to reveal singular behavior of tunneling current 
spectra in a high frequency region and to describe spectra 
peculiarities in a wide range of applied bias voltage. We 
found experimentally that changing of tunneling current 
noise spectra above the flat clean surface and above the 
impurity atom depends on the parameters of tunneling 
junction such as tip-sample separation or applied bias 
voltage 0. Our theoretical model is rather simple, more 
complicated models describing noise in semiconductors 
can be found in 

The investigations of the noise in two-level system was 
carried out in Authors studied current noise in a 
double-barrier resonant-tunneling structure due to dy- 
namic defects that switch states because of their interac- 
tion with a thermal bath. The time fluctuations of the 
resonant level result in low-frequency noise, the charac- 
teristics of which depend on the relative strengths of the 
electron escape rate and the defect's switching rate. If 
the number of defects is large, the noise is of the 1//" 
type. Shot noise in a mesoscopic quantum resistor is 
studied in @. Authors found correlation functions of 
all order, distribution function of the transmitted charge 
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and considered Pauli principle as the reason for the fluc- 
tuations. Current fluctuations in a mesoscopic conductor 
were analysed in Q. Authors derived a general expres- 
sion for the fluctuations in the cylindrical tunneling con- 
tact in the presence of a time dependent voltage. In 
noise of tunneling current at zero bias voltage was inves- 
tigated. It was demonstrated that at zero bias voltage 
the 1//" component of noise in the tunneling current 
vanishes and white noise becomes dominant. The 1//" 
dependence of the current noise in STM experiments on 
graphite in ambient conditions was investigated in 
Authors attributed this effect to fluctuations induced by 
adsorbates in tunneling junction area. In Q the fluctua- 
tions of tunneling barrier height have been investigated. 
The experiments have been performed under UHV condi- 
tions on graphite and gold samples using Ptir tips. From 
these measurements the authors have concluded that the 
intensity of barrier height fluctuations correspond to the 
intensity of tunneling current 1//" noise in the frequency 
range from 1 to 100 Hz. 

In our previous work we have found that for the low fre- 
quency region sudden switching on and off of additional 
Coulomb potential in tunneling junction area leads to 
typical power law dependence of tunneling current noise 
spectra at the threshold voltage In this article our 
aim is to study one of the possible microscopic origins 
of 1//" noise in tunneling contact in the non-resonance 
cases and to study high frequency peculiarities of tunnel- 
ing current spectra in a wide range of applied bias volt- 
age. We shall analise modiflcation of tunneling current 
noise spectrum by the Coulomb interaction of conduction 
electrons in the leads (metallic tip and surface) with non- 
equilibrium localized charges in tunneling contact. When 
electron tunnels from or to localized state the Coulomb 
potential is suddenly switched on or off and modifies tun- 
neling transfer amplitude. It will be shown that correc- 
tions to the tunneling vertexes caused by the Coulomb 
potential switching on and off result in nontrivial behav- 
ior of tunneling current noise spectrum in a wide range 
of applied bias voltage and should be taken in account. 
It was theoretically proved and observed experimentally 
that the same effects can lead to power law singularity in 
the current- voltage characteristics not only the threshold 
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voltage but also at the high frequency region of tunneling 
current spectra 



II. THE SUGGESTED MODEL AND MAIN 
RESULTS 

We shall analyse model of two localized states in tun- 
neling contact. In this case one of the localized states is 
formed by the impurity atom in semiconductor and the 
other one by the tip apex. 

When electron tunnels in or from localized state, the 
electron filling numbers of localized state rapidly change 
leading to appearance of localized state additional charge 
and sudden switching "on" and "off' Coulomb potential. 
Electrons in the leads feel this Coulomb potential. 

The model system (Fig. 1) can be described by hamil- 
tonian H: 

H = Hq + Htun + Hint 

^0 ^^i^p- eV)c+Cp + Skclck + ^ e^a+ai 

p k 2=1,2 

Htun = ^ TkiC^ai + ^ Tp^c+Gi + T ^ a^a2 + h.c. 

k,i p.i 

Hint = ^ Wiclcf^-aiol + V72C^Cfc'a2aJ 

Hq describes free electrons in the leads and in the local- 
ized states. Htun describes tunneling transitions between 
the leads through localized states. Hint corresponds to 
the processes of intraband scattering caused by Coulomb 
potentials Wi, W2 of localized states charges. 

Operators c^(cfc) and Cp (cp) correspond to electrons in 
the leads and operators af{ai) correspond to electrons in 
the localized states with energy e^. 




Current noise correlation function is determined as: 
{h/ef ■ S{t,t') ^< lL{t) ■ lL{t') >= 

= E n<4it'M'H{t)ci{t)> 

k,k 

where 

h{t) =Y,h.k-e= {J2c+it')adt')Tk, - h.c.) ■ | 

k k 

The current noise spectra is determined by Fourier trans- 
formation of ^(i, t'): S{uj) = J S{T)dT-e''^^. We shall use 
Keldysh diagram technique in our study of low frequency 
tunneling current noise spectra [ll|. 

Let's consider that STM tip is a metal with a cluster on 
the tip apex. In this case tip apex localized state energy 
coincides with tip Fermi level and consequently eV — £2- 

When we measure tunneling current spectrum above 
flat clean surface we deal with localized state connected 
with the tip. Tip apex localized state is usually formed by 
a cluster of several atoms or even by a single atom near- 
est to the surface. Tunneling transfer amplitude between 
the sample and the tip decreases exponentially with in- 
creasing of tip-sample separation. Thus in STM junction 
effective tunneling from or to the sample occurs through 
the atomic cluster nearest to the sample surface. So the 
tip apex cluster can be considered as an intermediate sys- 
tem for electron tunneling from the sample to the rest of 
the tip. Tip-sample separation is of order of interatomic 
distances. That's why the electron hopping from tip apex 
cluster to the bulk of the tip can be also treated as tun- 
neling process just in the same manner as it is done in 
tight-binding approach. 

In our case of weak interaction between localized states 
(T < 7fci,7fe2,7p2) possible variants for localized states 
energy levels position in tunneling contact are: 1. im- 
purity atom localized state energy level exceeds tip apex 
localized state energy level (ei > £2 = eV), 2. tip apex 
localized state energy level exceeds impurity atom local- 
ized state energy level (£2 = eV > ei). 

Expression which describes tunneling current noise 
spectra without Coulomb re-normalization can be cal- 
culated with the help of Keldysh diagram technic ■ It 
consists of three parts. 



FIG. 1: Schematic diagram of tunneling processes through 
states localized on impurity atom and on the STM tip apex. 



Soi{uj) and So2{(-u) are rather simple parts and they 
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FIG. 2: Lowest order diagrams contributing to tunneling cur- 
rent noise spectra for two localized states in tunneling contact, 
a) In the absence of Coulomb re-normalization of tunneling 
vertexes. b) In the presence of Coulomb re-normalization of 
tunneling vertexes. Tunneling vertexes are marked by the 
dot. Coulomb potential W2 is marked by the star. Solid lines 
correspond to electron Green functions. 



have the form: 



• ImGf,{uj + uj') ■ {n,{uj + uj') ~ I) ■ 

■ {ni{uj') - nk{uj')) + ni{uj') ■ {ni{uj + lo') 

- nu{uj + uj') + ■ [ du:'ImGf,{Lo') ■ 



■ ImGfi{uj + u') ■ {nk{uj + lu') - 1) • ni{uj') - 

- Uiiuj') ■ {ni{uj + uj') - 1) - rikiuj + lo') ■ nk{uj') + 

+ nk{uj') ■ {ni{LO + uj') + 'fki ■ J duj' ImGf^iuj + uj') ■ 

■ (ukiij')) ■ + uj') - 1) + ImGfiiuj') ■ 

■ (riiiuj')) ■ {nk{uj + uj') - 1) = So^ 



where i — 1 corresponds to the 6*01(0;) and i = 2 cor- 
responds to the 602 (w). Sosiuj) is not trivial, it exists 
only due to electron tunneling transitions from one lead 
to both localized states. Green functions shown on the 
graphs arc found in [l5l |. The contribution of Soi{uj) is 
given by graphs with i = j = 1 (Fig. 2a), So2{uj) is de- 
scribed by graph with i = j = 2 (Fig. 2a), 6*03(0;) is given 
by diagrams with i ^ j (Fig. 2a) . 



( ft/e)2 . So3(^) = 8 • 7fc, • 7fc, • J duj'ImGf^iuj') ■ 

■ ImG22{'^ + ■ {ni{uj') ■ {n2{uj + u}') - 1) + 
+ Hkiuj') ■ {Ji2{uj + uj') - 1) + {ni{uj') ■ 

■ {nk{uj + uj') - 1) - rikiuj') ■ {nk{uj + uj') - 1) + 

+ 8 • 7fc, • 7fe, • J duj' ImG^2{^') ■ ImG^^{uj + uj') ■ 

■ {n2{uj') ■ {ni{uj + lo') - 1) + 

+ Ukiuj') ■ {ni{uj + uj') - 1) + {n2{uj') ■ 

■ {nk{uj + uj') - 1) - nk{uj') ■ {nk{uj + uj') - 1) = 693 



Some typical low frequency tunneling current noise 
spectra for different values of dimensionless kinetic pa- 
rameters without Coulomb re normalization are shown 
on (Fig. 3a). 

It is clearly evident that when frequency aspire to zero 
tunneling current spectra aspire to constant value for dif- 
ferent dimensionless kinetic parameters. When impurity 
atom energy level exceeds tip apex localized state en- 
ergy level (ei > £2) we have the largest contribution to 
tunneling current spectra at the high frequency range of 
tunneling current spectra (w = £2 — £1). When tip apex 
localized state energy level exceeds impurity atom en- 
ergy level (£2 > £1) we have the largest contribution to 
tunneling current spectra at zero frequency. 

Now let us consider re-normalization of the tunnel- 
ing amplitude and vertex corrections to the tunneling 
current spectra caused by Coulomb interaction between 
both localized states and tunneling contact leads. Re- 
normalization gives us two types of diagrams contribut- 
ing to the final tunneling current noise spectra expres- 
sion. Ladder diagrams is the most simple type of di- 
agrams which gives logarithmic corrections to vertexes. 
But this is not the only relevant kind of graphs. We must 
consider one more type of graphs (parquet graphs) which 
also gives logarithmicaly large contribution to tunneling 
spectra. In parquet graphs a new type of "bubble" ap- 
pears instead of "dots" in ladder diagrams. In this situ- 
ation one should retain in the n-th order of perturbation 
expansion the most divergent terms (Fig. 2). 

It is necessary to re-normalize each vertex individually 
and re-normalize both vertexes jointly (Fig. 2b) p^ . 

The final expression for tunneling current noise spectra 
after Coulomb re-normalization of tunneling vertexes can 
be written as: 
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( h/ef-S{Lj) = So{Lj) + Soi{Lj) 
D2 
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{u + eV + Eif+Tl 



+ 5o2H-(( 
+ ( 
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{-uj + eV + Ei)^+rl 

£>2 



W2iy 



+ ( 



{-uj + eV + E2)^+Tl 

(( — + 

^^{io + eV + Ei)^+Tl' 



(( 



{iv + eV + E2)^ + rl 



+ ( 



{-w + eV + Ei)^+Tl 



(-W + eV + E2Y + r 



where, 



„ £l+£2 ^ V(£l-£2)^+4r2 
o = zL 

^'^ 2 2 

Ti ~ r2 ~ (7fc2 + 7p2 + 7/sl) 

where D- are the bandwidths of right and left leads, v- 
thc equilibrium density of states in the tunneling contact 
leads, W- Coulomb potential. 

In the case of strong interaction between localized 
states in tunneling contact we have localized states en- 
ergy levels splitting and as a result there is no singularity 
in tunneling current spectra at the low frequency region. 
In our situation of weak interaction between localized 
states (T < 7fci,7fe2,7p2) tunneUng transfer amplitude T 
plays an important role in kinetic processes in the tunnel- 
ing junction but weakly influences on the energy spectra. 

First of all we consider the situation when both lo- 
calized states acquire positive charge. Fig. 3b demon- 
strate low frequency tunneling current noise spectra for 
typical values of dimensionless kinetic parameters. We 
can sec that re-normalization of tunneling matrix ele- 
ment by switched "on" and "off' Coulomb potential of 
charged impurities leads to typical power law singularity 
in low frequency part of tunneling current noise spectra 
and to the peak in the particular high frequency region, 
caused by the singularity on the frequency a; = £2 — £1. 
Peak amplitude depends on kinetic parameters and on 
the charged localized states Coulomb potentials. 

At the fixed composition of kinetic parameters max- 
imum peak amplitude corresponds to the case when 
charged impurity atom localized state Coulomb poten- 
tial W\ exceeds Coulomb potential of tip apex localized 
state W^- In the case of W2 > W\ (Coulomb potential 



of charged impurity atom is larger than Coulomb po- 
tential of charged tip apex localized state) peak ampli- 
tude slightly changes by comparison with the case when 
Coulomb potential of charged impurity atom is similar 
to Coulomb potential of charged tip apex localized state 
Wi - 14^2- 

Tunneling current noise spectra in double logarithmic 
scale demonstrate frequency regions where every part of 
final expression which include different power law expo- 
nent, approximate noise spectra in the best way Fig. 3c. 

Let's analyse tunneling current spectra shown on 
Fig. 3c-e. In the case of two interacting positively charged 
localized states the tunneling current spectra at low fre- 
quency and in the region of the high frequency singu- 
larity is always determined by the term which produces 
the most strong of logarithmic singularity, determined by 
the sum of localized states Coulomb potentials at any val- 
ues of dimensionless tunneling rates of tunneling contact. 
(Fig. 3c). 

If one of the Coulomb potentials strongly exceeds 
the other one this potential determine tunneling current 
noise spectra with increasing of frequency. If Coulomb 
potential of charged impurity atom is similar to Coulomb 
potential of charged tip apex localized state Wi ~ W2, 
tunneling current noise spectra with increasing of fre- 
quency is determined by Coulomb potential of charged 
tip apex state W2 (Fig. 3c). 

Let us consider tunneling current noise spectra in 
double logarithmic scale in the situation when localized 
states acquire charges of different signs Fig. 3d. 

Tunneling current noise spectra in double logarithmic 
scale (Fig. 3d) make it clear that when localized states 
acquire charges of opposite signs tunneling current noise 
spectra in the low frequency region and in the region of 
the high frequency singularity (cj = £2 — £1) are always 
approximated by the term depending on maximum pos- 
itive value of impurity atom Coulomb potential Wi or 
tip apex localized state Coulomb potential W2- Tunnel- 
ing current noise spectra in double logarithmic scale in 
the situation when both localized states acquire negative 
charges (Fig. 3e) make it clearly evident that singular 
effects become negligible. 

Now let's describe interaction effects in tunneling con- 
tact when both localized states are formed by impurity 
atoms in the surface. In this case we can say about clus- 
ter in the surface eV 7^ £2 7^ £1- 

Possible variants for localized states energy levels po- 
sition in tunneling contact are: 

1. applied bias voltage exceeds energy levels of both 
localized states in tunneling contact {eV > £2 > ei, eV > 
£1 > £2); 2. applied bias voltage exceeds energy level 
of one of the localized states in tunneling contact (£2 > 
eV > £1, si > eV > £2); 3. energy levels of both localized 
states in tunneling contact exceed applied bias voltage 
(£2 > £1 > eV, £1 > £2 > eV) 

Typical low frequency tunneling current noise spectra 
for different values of dimensionless kinetic parameters 
without Coulomb re-normalization are shown on Fig. 4a. 
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FIG. 3: Typical tunneling current noise spectra for different values of dimensionless kinetic parameters for two localized states 
in tunneling contact {eV = £2 7^ ei )• a) In the absence of Coulomb re-normalization of tunneling vertexes. b) In the presence 
of Coulomb re-normalization of tunneling vertexes. c) In the presence of Coulomb re-normalization of tunneling vertexes in 
double logarithmic scale when both localized states acquire positive charge, d) In the presence of Coulomb re-normalization 
of tunneling vertexes in double logarithmic scale when localized states acquire charges of different signs, e). In the presence of 
Coulomb re-normalization of tunneling vertexes when both localized states acquire negative charge. 



It is clearly evident that when frequency aspire to 
zero tunneling current spectra aspire to constant value 
for different dimensionless kinetic parameters. When 
energy levels of both localized states in tunneling con- 
tact exceeds applied bias voltage (£2 > ei > eV or 
El > 62 > eV) we have the largest contribution to tunnel- 
ing current spectra in the high frequency region (Fig. 4a) . 
In the other cases we have the largest contribution to 
tunneling current spectra at zero frequency (Fig. 4a). 

Now let us consider re-normalization of the tunnel- 
ing amplitude and vertex corrections to the tunneling 
current spectra caused by Coulomb interaction between 
both localized states and tunneling contact leads. Fig. 4b 
demonstrate tunneling current noise spectra for typical 
values of dimensionless kinetic parameters. It is clearly 
evident that when frequency aspire to zero tunneling cur- 
rent spectra aspire to constant value and we have no 
power law singularity in a low frequency part of tunneling 
current spectra. 

We can see that re-normalization of tunneling matrix 
element by switched "on" and "off' Coulomb potential 
of charged impurities leads to the singular peaks in high 
frequency regions of tunneling current spectra, caused 



by singularities on the frequencies to = eV — £\ and uj = 
eV-e2. 

Let's start from the situation when both localized 
states acquire positive charge. Tunneling current noise 
spectra in double logarithmic scale demonstrate fre- 
quency regions where every part of final expression which 
include different power law exponent, approximate noise 
spectra in the best way Fig. 4c. 

Let's analyze tunneling current spectra shown on 
Fig. 4c-e. In the case of two interacting positively charged 
localized states the tunneling current spectra in the re- 
gions of the singular peaks is determined by the term 
which produces the most strong of logarithmic singular- 
ity, determined by the sum of localized states Coulomb 
potentials at any values of dimensionless tunneling rates 
of tunneling contact. (Fig. 4c). 

If one of the Coulomb potentials strongly exceeds 
the other one this potential determine tunneling current 
noise spectra except frequency regions in the vicinity of 
singular peaks . If Coulomb potentials have the similar 
values Wi W2, tunneling current noise spectra is de- 
termined by Coulomb potential of charged tip apex state 
W2 (Fig. 4c). 

Let us consider tunneling current noise spectra in 
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FIG. 4: Typical tunneling current noise spectra for different values of dimensionless kinetic parameters for two localized states 
in tunneling contact {eV 7^ £2 7^ ei)- a) In the absence of Coulomb re-normalization of tunneling vertexes. b) In the presence 
of Coulomb re-normalization of tunneling vertexes. c) In the presence of Coulomb re-normalization of tunneling vertexes in 
double logarithmic scale when both localized states acquire positive charge. d)In the presence of Coulomb re-normalization 
of tunneling vertexes in double logarithmic scale when localized states acquire charges of different signs. e)In the presence of 
Coulomb re-normalization of tunneling vertexes in double logarithmic scale when both localized states acquire negative charge. 



double logarithmic scale in the situation when localized 
states acquire charges of different signs (Fig. 4d). 

Tunneling current noise spectra in double logarithmic 
scale (Fig. 4d) make it clear that when localized states 
acquire charges of opposite signs tunneling current noise 
spectra is always approximated by the term depending 
on maximum positive value of impurity atom Coulomb 
potential Wi or tip apex localized state Coulomb poten- 
tial W2. 

Tunneling current noise spectra in double logarithmic 
scale in the situation when both localized states acquire 
negative charges (Fig. 4e) make it clearly evident that 
singular effects after Coulomb re-normalization become 
negligible. Simple estimation gives possibility to analyze 
the validity of obtained results. For typical composition 
of tunneling junction parameters in the low frequency re- 
gion power spectrum of tunneling current corresponds to 
experimental results. Power spectrum on zero frequency 
has the form: 

^(0) « (7e// ■ e/nf ■ {Dh.ff.r'^ ■ (1/ A uj) 

For typical 7^//, 7em « 10"", D ^ 10, W ^ 0,5 

5(0) « 10-18 AVffz g. 



III. CONCLUSION 

The microscopic theoretical approach describing tun- 
neling current noise spectra in a wide range of applied 
bias voltage taking in account many-particle interaction 
was proposed. When electron tunnels to or from local- 
ized state the charge of localized state rapidly changes. 
This results in sudden switching on and off of additional 
Coulomb potential in tunneling junction area, and leads 
to singular behavior of tunneling current spectra in a par- 
ticular frequency range determined by the parameters of 
tunneling junction such as applied bias voltage or the 
localized states energy levels deposition. 

In the case of two localized states in tunneling junction 
when energy level of one of the localized states is con- 
nected with the tip apex and is not equal to energy level 
of the impurity atom localized state on the surface we can 
see typical power law dependence for the low frequency 
part of tunneling current spectra and singular peak in 
the particular high frequency region of tunneling current 
spectra determined by the localized states energy levels 
deposition in the tunneling junction (eV — e-i ^ ei). 
In the non-resonance case of the cluster on the surface 
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(eV =/= S2 ^ Si) wc have found two singular peaks in the 
high frequeney region of tunneling current noise spectra. 
So our results demonstrate that changing of the applied 
bias voltage leads to the tuning of tunneling current noise 
spectra . 

This effect can be qualitatively understood by the fol- 
lowing way: apart from the direct tunneling from the 
localized state to the state with momentum k in the lead 
of tunneling contact described by the amplitude Tk, the 
electron can first tunnel into any other empty state k' 
in the lead and then scatter to the state k by Coulomb 



potential. So the appearance of even a weak Coulomb 
potential significantly modifies the many particle wave 
function of the electron gas in the leads. 
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